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Abstract 



We present an approximate model of Wheeler-Feynman electrodynamics for which unique- 
ness of solutions is proved. It is simple enough to be instructive but close enough to Wheeler- 
Feynman electrodynamics such that we can discuss its natural type of initial data, constants of 
motion and stable orbits with regard to Wheeler-Feynman electrodynamics. 
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1 Introduction 

Already in the early stages of classical electrodynamics it was observed that the coupled equations 
of motion of Maxwell and Lorentz for point-charges are ill-defined. Since then many attempts have 
been made to cure this problem among which the most famous one is Dirac's mass renormalization 
[Dir38]. However, none of the known cures has yet led to a mathematically well-defined and phys- 
ically sensible theory of relativistic interaction of point-charges. Rather than being a cure of the old 
theory, Wheeler-Feynman electrodynamics (WF) is a new formulation of classical electrodynamics 
and is maybe the most promising candidate for a divergence-free theory of electrodynamics about 
point-charges that is capable of describing the phenomenon of radiation damping [WF45, WF49]. 
For an overview on the mathematical and physical difficulties of classical electrodynamics and the 
role of WF we refer the interested reader also to the introductory sections of [BD01] and [BDD10]. 
In contrast to text-book electrodynamics, which introduces fields as well as charges, WF is a theory 
only about charges - fields only occur as mathematical entities and not as dynamical degrees of 
freedom. Let N be the number of charges, q i t e K 3 the position of the z-th charge at time t, nit > 
its mass, and 




dt 



its relativistic momentum and velocity. The fundamental equations of WF take the form 



dt \Pi,t 




1 < i < N. 



(1) 
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We have chosen units such that the speed of light equals one. The force term F,[q 7 ] is a functional 
of the trajectory t h-» q ; and can be expressed by 

F f [q ; ](x, p) := e } £ (Ef[q,-](x) + v(p) A Bf [q ; ](x)) 



where (E+[q 7 ],B+[q 7 ]) and (E~[q 7 ],B^"[q 7 ]) denote the advanced and retarded Lienard-Wiechert 
fields - in our notation E represents the electric and B the magnetic component and A is the outer 
product. The Lienard-Wiechert fields are special solutions to the Maxwell equations corresponding 
to a prescribed point-charge trajectory 1 1-» q 7 - jf ; see e.g. [BDD10]. Their explicit form is 



[ (n 7 ±v ; )(l-v2) 

Ef[q,](x) := - 2 1 + 

L||x-q,|| (l±n r v y ) 3 

Bf[q 7 ](x) := +n /± AEf[qj](x). 
Here, we have used the abbreviations 

x - q/.r dq u 



n j A ( n ; A a ;j 
|x-qJ|(l ±n r v 7 -) 3 



(2) 



x - q/ri 



v ; = 



dt 



d 2 <lj,t 



f=f± 
j 



dt 2 



(3) 



and the delayed times t + . and t- are implicitly defined as solutions to 



tj(t,x) :-- 



t + 



X - Hj,tj(t,x) 



(4) 



Given the space-time point (t, x) the delayed times ft and fT are determined by the intersection 
points of the trajectory t h-> qj t with the future and past light-cone of (t, x), respectively. As long 
as k/n < 1 both times ?t and ?t are well-defined. As it becomes apparent from (4) the WF 
equations (1) involve terms with time-like advanced as well as retarded arguments which makes 
them mathematically very hard to handle: 

• The question of existence of solutions is completely open with the exception of the fol- 
lowing special cases: Schild found explicit solutions formed by two attracting charges that 
revolve around each other on stable orbits [Sch63]. Driver proved existence and uniqueness 
of solutions for two repelling charges constrained to the straight line whenever initially the 
relative velocity is zero and the spatial separation is sufficiently large [Dri79]. Furthermore, 
Bauer proved existence of solutions in the case of two repelling charges constrained to the 
straight line [Bau97], and a more recent result ensures the existence of solutions on finite but 
arbitrarily large time intervals for /V arbitrary charges in three spatial dimension [BDD10]. 

• It is not even clear what can be considered to be sensible initial data to uniquely identify 
WF solutions. While Driver's uniqueness result suggests the specification of initial positions 
and momenta of the charges, Bauer's work points to the use of asymptotic positions and 
velocities to distinguish scattering solutions, and one sentence below Figure 3 in [WF49] 
hints to a certain configuration of whole trajectory strips of the charges as initial conditions. 
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Furthermore, one important issue of WF is the justification of the so called absorber condition and 
the derivation of the irreversible behavior which we experience in radiation phenomena [WF45]. 
The link between WF and experience must be based on a notion of typical trajectories, i.e. on a 
measure of typicality for WF dynamics. Such a measure is unknown and at the moment out of 
reach. A generalization to many particles of the approximate model we consider next provides an 
excellent case study for introducing a notion of typicality for this kind of non-markovian dynamics. 
Note that uniqueness as well as conservation of energy are very likely to be important for defining 
a measure of typicality in the sense of Boltzmann. This is work in progress. 

Our intention behind this work is to provide a pedagogical introduction to the mathematical 
structures arising from delay differential equations of the WF type by discussing initial data and 
uniqueness, constants of motion and stable orbits in a non-trivial approximate model of WF, which 
was introduced in [Dec 10], and for which many mathematical questions can be answered satisfac- 
torily. Compared to (1) we make the following simplifications: 

• We consider only two charges, i.e. N = 2. 

• As fundamental equations of this approximate model we take the form (1) where we replace 
the Lienard-Wiechert fields (Ef[q y ](x), Bf [q y ](x)) with 

E^tqyKx) := ej '- 3, Bf [q,] = 0, (5) 

X - q;,f±(f, X ) 

i.e. the Coulomb fields at the respective delayed times in the future and the past. 
The resulting equations are 

d /q, A / v(p i>( ) \ 

5 ip,,) = U[F(q M -,^„ ) ) + F(, i ,, -qvpJjQ- *'«<••* >* <6) 
where the force field is given by 

F:R 3 \{0}^R 3 \{0}, x ^F(x):=^. (7) 

l|x|| 

We emphasize that the delay function (4) is the same as in WF and that the simplified fields (5) 
are the longitudinal modes of the Lienard-Wiechert fields (Ef [q y ](x), Bf [q y ](x)), i.e. 

V • (Ef [q 7 ](x) - E^tq^-Kx)) = 0, V A E^Kx) = 0. 

Furthermore, for small velocities and accelerations of the j-th charge one finds Ef [q ; ](x) « E['' ± [q ; ](x). 
In this sense one can regard this simplified model as an physically interesting approximation to WF. 
Note also that, in contrast to Ef [q 7 ], the simplified field Ef'^qy] does not depend on the accelera- 
tion of the j'-th charge; compare (2). This fact makes it easier to control smoothness of solutions, 
however, is not the key difference that allows us to prove uniqueness of solutions for the approxi- 
mate model. 



This paper is organized as follows: 
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1. In Section 2 we discuss natural initial data for the approximate model (Definition 2.1), and 
show how to construct solutions uniquely depending on given initial data (Theorem 2.3). A 
byproduct is the observation that in general the specification of initial positions and momenta 
of the two charges is not sufficient to ensure uniqueness (Corollary 2.6). 

2. In Section 3 we discuss the structure of constants of motion by reference to the energy 
conservation in the approximate model. The energy functional is defined in Definition 3.1 
and the energy conservation is shown in Theorem 3.2. 

3. In Section 4 we identify stable orbits (Definition 4.1) which solve (6) as proven in Theo- 
rem 4.2. 

4. We conclude in Section 5 by putting these results in perspective to WF. 



• If not otherwise specified we use the convention that i, j e {1,2} and j ^ i. 

• Any derivative at a boundary of an interval is understood as the left-hand or right-hand side 
derivative. 

• Vectors in x, y e R 3 are denoted by bold letters, the inner product by x • y, the outer product 
by x A y, and the euclidean norm by ||x||. 

• V, V-, and VA denote the gradient, divergence and curl w.r.t. x, respectively. 

• Overset dots denote derivatives with respect to time, i.e. q, , = ^ and = ^p*. 

• The future and past light-cone of the space-time point (t, x) is defined as the set 



for + and -, respectively. 

2 Construction of Solutions 

For this section we regard the indices (i, j) fixed to either (1, 2) or (2, 1). We consider the following 
initial data (see Figure la): 



Definition 2.1. We call a pair (q: , j _ of smooth position and momentum trajectory stripes 



Notation. 



L ± (t, x) : = {(s, y)\(t- s) 2 - (x - y) 2 = 0, t > ±s} 




(0) . 



R D [tf\ ^] 



R- 



3 



^q 



(0) 



(0) . 



R- 



3 



I €{1,2} 



that fulfill the following conditions initial data for equation (6): 



(i) For times t e [tf\ <J 1} ] n [tf\ tf] it holds 
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(a) The thick black stripes represent the ini- (b) Nature of the interaction, 

tial data. 

Figure 1 



(ii) For all t e [tf\ tf } ] it holds 

a n (0) 

§<C = v(p»») ^ —M==, (.,1.2,. (8) 

V m * + \At II 

(iii) The times tf\ tf \ t ( p, relate to each other according to 

,f> = tt (<f,q^), .f^fU™), ^ ('«''• <»)■ W 



(iv) At time t = f| 0) and for all integers n > the trajectories obey 

<F_ /q, A = cT_ t v( Pi - 1 ) 

<fc" \p, J ^ U<v [F (q,-, f - q./.ru,,..,) + F (q M - q ; ., „.„..,)] 

(v) At time t = f^and for all integers n > the trajectories obey 

dt" \Pj,t) dt" [ejei [F (q p - q^ ((>M ) + F (q p - q,-, f r (f , q ,,))]/ " 
Remark 2.2. Note that (8) requires ||q !jf || < 1 for i 6 {1,2} and 



Pi, 



Furthermore, such initial data can be constructed as follows: 



(10) 



(11) 
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1. Choose times z^ 0) < tj~\ an arbitrary smooth trajectory strip q^ 0) : [tj^jj] — > K 3 with 
q^.j < 1, one space-point (f: , q. ( 0) ) on the intersection of the forward light-cone of 



q L) with the backward light-cone of (r. , q ,,,), and one space-point (f. , q L) some- 

j.t. "~ } /,1 1 ' IX 

where on the forward light-cone of (/p, q (0) (1) ) and inside the forward light-cone of (tf\ q (0 ( 0) ). 



(0) ._ "'jC 



2. Define p^/ := p= for f 6 [f- ,r. ] and, using (6), compute all derivatives of p i? at the 

* _ /o) „„ , . _ ,(1) 



times t = ?; 0) and t = t\ 



3. Choose a smooth trajectory strip q^ 0) : [t- ,t.] 



I 3 through the space-time points (tf\ q (0 L) 

1 i.t. 



and (tf\q (0) w ) such that 



m JO) 

q^ < 1 and that t i-> := , "' '•' smoothly connects to the 

VHM 2 



derivatives computed in step 2. 
Provided such initial data we shall prove our first result: 
Theorem 2.3. Given the initial data (qf\ pP) _j 2 there exist two smooth maps 



2 D, 



x 



(12) 



W (q,-. f , Pi,;) = (q-?, pS?)/or all t e [t?>, t) 1J ] and i e {1,2}. 



(0) ,(l) n 



(«J * h> (q,-, f , p,-, f ) so/ves (6J on D ; := {l*{Tf* t q^), ^(r™ x , q^))/or f e {1,2}, / * i. 



(nij For j = 1,2 let D, c E, be an interval such that {f, ,t. ] c D, an J Ze? Z), — > 
? i — ► (q i f , p i f ) oe a smooth map such that t h-> (q,- ,, p,- ,) solves (6) for t e D,-. r/?en: 

(q /if ,Pi >? ) = (q,>p ; -, f ) W e D, n D„ z e {1,2} 

(q^,Pi, f ) = (q^P^W 6 [tf\tf\i e {1,2}. 

Given appropriate constants d > an J < v < 1, the times -oo < r™ n < zj 0) < f: < r ; max < oo 
are defined such that \T™ m , r max ] is the largest interval containing [tf \ t\ ] wif/i f/ie property: 



I 3 x R 3 , 



q;,f - q/ff (?,q,,,) 



><*, e [r ( min ,r, max ] 



(13) 



Their value is determined during the construction of (12) in the proof. 

Remark 2.4. Our focus lies on the uniqueness assertion (iii) of Theorem 2.3. We do not attempt 
to give a priori bounds on T™ n , T™ ax , i e {1,2}, whose values are determined during the dynamics 
by condition (13). This condition is needed to prevent two types of singularities that can occur: 
First, the approach of the speed of light, and second, collision or infinitesimal approach of charges. 
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These singularities can also be present in WF which can be seen directly from the form of the fields 
(2). While the second one is familiar since it is of the same type as seen in the Af-body problem of 
Newtonian gravitation [SMSM71], the first one is very specific to WF-type delay problems. Such 
singularities are due to the nature of the delay times fj and t~ defined in (4) which tend to plus or 
minus infinity if the j'-th charge approaches the speed of light in the future or the past, respectively; 
the origin of this singularity can be seen best in (18). Because of angular momentum conservation 
it is however expected that fox N = 2 charges one always finds T min = -oo and T max = oo, which is 
at least true for the solutions given in Section 4. A treatment of the Af-body problem will require a 
notion of typicality of solutions. 

The key ingredient of the proof, which can be checked by direct computation, is the following: 
Lemma 2.5. The following statements are true: 

(i) The map F defined in (7) is bijective and its inverse is given by 



I 



\{0} 



\{0}, 



y i-> i(y) - 



|3/2 • 



(ii) Let (q| 0) ,P- 0) )._ be given initial data. For any integer n > equations (10) and (11) are 
equivalent to 

^ (q* - q/.r „„..,) = ^1 \^Vut ~ F (q,v - q,ru M ^ for t = tf\ 



and 



d n d n I 1 \ 

^ {qj, t ~ q^J) = ^1 y— Pjt ~ F (q ;> - q^^) j for t = tf, 



respectively. 

Lemma 2.5 ensures that for example in situations as depicted in Figure lb we can compute 
from F(q, ;f - q ; , ; +), which is determined by the initial data and (6), the space-time point (ft, qp+). 
This is the key ingredient in our construction: 



Proof of Theorem 2.3. As a first step, we construct a smooth extension of q^ 0) beyond time t y } l> . Let 
us introduce the short-hand notation 

tj = tj(t,qij). 

In general, any solution t h-> (q i>t , p,v)i=i,2 to (6) has to fulfill 



Pi t = eiej [F (q u - q j>t j) + F (q u - q p -)] . 
With the help of Lemma 2.5(i) we can bring this equation into the form 

qj,t+ = q*, t - I (^P^ - F (q/,r - q^Jftq,,))) 
for times t e [r ; min , r ; max ] because then 

l(^PM-F(q*-q* 7 )) 



(14) 



(15) 



U,t Hj,tj 



> d 
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is guaranteed, and hence, the right-hand side of (15) is well-defined. We now make use of (15) to 
compute 



according to 



t] = t + 



q (1) + 



t| 1 A (0) v(„(P) „(0) \\ 
_(0) t/ 1 A (0) W /_(0) _(0) \| 



(16) 
(17) 



for all t e (/f'./J 1 ^ n (r ; min , r ; max ). Note that due to (9) the right-hand side of (16) and (17) is 
well-defined. We define 

+(2) _ • „ Jy.+ i T'max I 

f. -mm |^Ua),^ ]■ 
Since q| 0) is smooth on , t\ ] also ft depends smoothly on? 6 ,^ ], and in consequence, 



qf is smooth on (t\ l> , f?']. Furthermore, a direct computation gives 



AT) 



t H> 



1 + n, + • qg 
l+^ + -q^' 



(0) 


-Q (1) 


i,t 




(0) 
i,f 


-Q (1) 





(18) 



so that, using the notation 



d dt d 
dfj dfjdf 



we may then compute 



lim — q (1) 

. 5V <» ds"^ s 



d n d n 
lim -r— q (1 ?+ = lim — — 



„(0) t/ 1 A (0) W /„(0) u 



(19) 



for every integer n > 0. Lemma 2.5(ii) ensures that 



(0) 



and hence, 



t i-> 



q^ far fe [if, if] 
qf for I 6 (if, if] 



is a smooth map on [if, if). Furthermore 



Vie (if, if]. 



(20) 



is well-defined by (13). 
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In the second step, we use the analogous construction to extend q^ 0) smoothly beyond time tf: 
We define 



(i). 



by 



tj = t + 



l j,t 



I 



F q 



q (0) 



for all t e (tf, tf] n [7/™ n , 7/ ; max ] and furthermore 



♦(2) _ . ■ „ i Tmax I 

t. - rnin^,. \ t= pJi j. 



As in the first step one finds that 



\qf! for ? e[f,f] 
[qg> for ^ C^,^] 

is smooth for ? e [^ 0) , ?| 2) ]. Finally, due to (13) we can define 



(0) Jlh 



1 1 — ► 



p; 



(i) .. 



flu 



m :, 



In consequence, the maps 



for i = 1,2 are smooth, and by virtue of definitions (20),(23) and (17), (22) and (10),(11) of (2.5) 
they fulfill 



,(0) 



(21) 
(22) 



(23) 



for?e[?f ) ,?r(^ ) ,q. t c 2) )]and 



P;,f 



v(Py) 

^[F( qi , t -q^) + F(q i)( -q 7;(7 )] 



V (P;>) 

eyei [F(q^-q ; - ir )+F(q y , t -q Mr )] 



(24) 



(25) 



for/G^.fT^.q.^)]. 

This construction can be repeated where in the kih. step one constructs the extension 



( (k) (4+1)-. . , (*) (4), 



For each step one finds 

tf +1) > min {t (k) + d, r, max ) , i G {1,2}. 
In consequence, only finite repetitions of this construction are needed to compute 



|(q!?,pi?) for te[tf\t^] 
l(q*.pff) forfeit?,^ 



[tf\ r- x ] : 1 i > (q (> p,-,) := < ^''g' " ~ ^' (26) 



which fulfills (24) for f e [? f,t7(Tf^,q iT ^)] and (25) for f e [t) VTCTP^^, ^)]. The same 

l I J J* j J J I ' i 

construction can be carried out into the past which results in smooth maps 

[7f n , ?T X ] -> R 3 x R 3 , t» (q ;> p u ), i e {1, 2}. 

From this construction we infer the claims of Theorem 2.3: Claim (i) follows from definition 
(26). Furthermore, due to Lemma 2.5, (16)-(17) and (21)-(22) the map t h-> (q^, p (j ,) fulfills (6) for 
times 

t e [tt(Tf a , q jJr ), t-{T™\ q u: .~)\ 

for i 6 {1, 2} and j ^ i and therefore claim (ii). Finally, Lemma 2.5 guarantees that this constructed 
solution is unique which proves claim (iii) and concludes the proof. □ 

As a byproduct we observe that specification of initial positions and momenta of the charges 
as suggested by the work [Dri79] does not always ensure uniqueness: 

Corollary 2.6. Let (qf\ P- 0) )._j be initial data such that there is a t* e (tf\ f ; a) ) n (Sp, t ( p), let 

Dj — > R 3 xR 3 ,f i-> (q i>t ,pi tt )fori e {1,2} be the corresponding solution to (6), and let qf } ,p| 0) 6 R 3 
for i = 1,2 be defined as 

(q^pf^ta-V'P^L' »'e{l,2}. 

77zere are uncountably many other solutions D { — > R 3 x R 3 , 1 1-» (q,-, f ,p,-, f )/or z = 1, 2 to ((5) which 
fulfill 

(q%Pv)L =(q! 0) ,pS 0) ), i€{l,2}, (27) 

pm? no? 

(q ! >P,r) = (q 1 >p 1 -, f ) VfeAnA, ie{l,2}. (28) 

Prao/ Choose 5 e R and 5 > such that (j - 5, s + 5) c (zf \ zf } ) n (/ 0) , ^) and f* £ (j - tf, $ + 5). 
Furthermore, for A > and z = 1, 2 let d^,- : R — » R 3 , ? i-» d^Cf) be a smooth function such that 

supp d AJ = [s - 6, s + 5], sup |da,,(OI = A. 

We define 

choose the parameter A > such that ||q, ,|| < 1, and define 

Vi - 
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\ 

u 

t + (t.) 


A 


1^3, t 
t + (t,) 

t~(ti) 














\ \ 















Figure 2: Given t t , tj the thick lines denote the data required to define the energy of the system. 

The maps (qj 0) , p| 0) ),=i,2 are initial data according to Definition 2.1 which fulfill (27). However, 
according to Theorem 2.3 the solution D, — > (r 3 x R 3 ),? i-> (q,> p,-,r) for i = 1,2 to (6) corre- 
sponding to (q| 0) , P; )i=i,2 does not fulfill (28). Note that there are uncountably many choices, e.g. 
in s, 6, A and d^,-, to define other (qj ,p) )/=i,2 such that (27) holds. None of the corresponding 
solutions however fulfill (28). □ 



3 Constants of Motion 

In the following we define an energy functional for the approximate model from which the general 
structure of constants of motion will become apparent. Throughout this section we consider a 
solution 

R □ D t -> R 3 x R 3 , (q,- „ Pi ,), t e {1, 2}, 

to (6); see Theorem 2.3. 

Definition 3.1. We define a map // : D! x D 2 — > R + by 

2 . j 2 1 

where for the delay functions we have used the short-hand notation 

tf(t) = tf(t,q u ), ie{l,2}, j * i. 

We refer to H as the energy functional of the system. See Figure 2 for an example of which data is 
needed to define this functional. 
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Theorem 3.2. For all t\ e D\ and t 2 e D 2 the equality 



H(t u t 2 ) = H(0,0) 

holds true. 

Proof. Since we assume a smooth solution t i-» (q, ,, p,,r),=i,2 to (6) it is straight-forward to verify 
the claim directly by differentiating with respect to ti and t 2 . Here, however, we want to provide 
a general idea how to find constants of motion for a WF-type delay differential equations, and 
therefore give a more instructive proof: 

We start with the sum of the kinetic energy difference between times and f, of the two particles 
i = 1,2, that is 

V dsp ks -\($ ks ) (29) 
ti J ° 

and make use of the equations (6) to express this entity in terms of 

(29) = 2_ j ei2_ J e j2-j dfF (*.*-^w)-*> (30) 

i'=l j*i ± J° 

It is convenient to split (30) into the following summands: 

1 C ! ' 

< 30) = oZ e <Z e jZ dsF(q Us -q ms) )-q Us (31) 

i ^ r tl l dt% s y 



1 r'' / a? (^)\ 

+ 2 Z e; Z e ; Z I J * F(q <> " q ;-'> ,) • fc* - -tt~ (32) 

i'=l ± ^° ' ' 

i r tf 

+ 2 Zj ei Zj Zj ^ F(q '- 4 ' ~ q ^ ±w) ' q/ ° w_ 7T- (33) 

i'=l Mi ± ^° 



(32) = 4 Z e <- Z e ; Z " — 1 — + c (34) 



j* 

The integrand in (32) is an exact differential so that 

2 

2 '— 

i=i & ± ||q/, fj - q/^fe) 

where C e R is a constant. Next, we exploit the symmetries of the force field and the delay function 

F(x) = -F(x), V x e R 3 \ {0}, (35) 

t = tj(tj(t)), ie{l,2}, (36) 
Now (36) allows to rewrite (33) by substitution of the integration variable according to 

(33) = n Z e ' Z e i Z I Js F (q^« - q/>) • q*> (37) 

z Vt j7i ± J 'J(0) 
12 



Furthermore, we apply (35) and after that relabel the indices i ±^ j and + ±^ + to get 

2 _ _ rite) 



1 \ i ptjw) 

(37) = ~ i Z e/ Z e J Z J * F(q ^' ~ q ^f ■ q ^ 

= _ oZ^'Z^'Z '/.vF(qa -q/./:M i )-(lo 

i 2 r° 
= _ oZ^Z^'Z ^ F ( q '>- q 7^w)- q '> 

_ o Z e ' Z e ^ Z ds F(q '- 1 ' ~ q ^> i) • q '> 

2 tl #i V J « 

~ \ Z g ' Z e > Z f F(q ^ ~ q ^'; ±($)) ■ 



(38) 
(39) 
(40) 



Noting that term (38) is just another constant and term (39) cancels the first term on the right-hand 
side of (31), the kinetic energy can be written as 



2 2 



(29) = ~ e ' Z e i Z ij — - \ Z e < Z e ; Z f Js F(q <> " q ;4 w) ■ q ^ + c 

i=l ± ||q,-r, - qy,f*(fi) t=l #i ± J '< 

which proves the claim. □ 

4 Stable Orbits 

By the symmetry of the advanced and retarded delay it is straight-forward to construct stable orbits 
(compare [Sch63]): 

Definition 4.1. For masses mi,m 2 > 0, charges ei,e 2 e R, radii ri,r 2 > 0, and angular velocity 
wel\|0( such that 



mi y(cur 2 )r 2 „ AT/ Jf ^ +r 

e!e 2 <0, — = — --, 0<wAr<-, AT := — 

m 2 y((on)ri 2 ^ 1 + -^nnyiujr^na) 2 

where 

1 

y(v) : 



we define the particle trajectories 



cosojt\ cos cot\ .... 

fi->qi ( :=ri . , t h-> q 2r := -r 2 . (41) 
M ' isinwf/' z \sinoW v y 



which we call Schild solutions; see Figure 3. 
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Figure 3: Stable orbit of two charges revolving around the same center with radii r ( , rj. AT denotes 
the modulus of the delay. 



Theorem 4.2. Given m,\m 2 , r\ > 0, e\ e R there are r2 > and co £ R such that the Schild 
solutions (41) obey (6). 

Proof. We start by rewriting (6) as a second order equation, i.e. 

j t {myiAiMi^ = e i e j [ F (fe - q^t&q,,,)) + F (qm - ( i / ? ./ <i 01 > z e ! 1 » 2} ' ' J- 

If the charges stay on the circular orbits (41) then velocities and accelerations are constant and 
fulfill 



\q iit \\ = co r t . 



Furthermore, the kinematics dictate that the net force acting upon the particles must be centripetal 
w.r.t. the origin and equal, i.e. 



(miy(qi, f )qi >f ) 



= m\y{u)r\)oj L r\ = m2y(cor2)to 2 r 2 



which implies 



m x y(cor 2 ) r 2 



m 2 y{tor{)ri 

Next, we compute the delay function AT > which, according to (4), was defined by 
AT 2 := (tj(t) - if = ||q I>r - q^JI =r 2 l +r 2 2 + 2r x r 2 cos (coAT) . 



(42) 



(43) 



Due to the symmetry in the advanced and retarded delay (see Figure 3) the net force is centripetal 
w.r.t. the origin for 

e x e 2 < 



and its modulus equals 

epj [F(q,- f - q ;, ,/,,..,) + v(q itl - q;, t T (f)qi , () )]| = 

14 



e x e 2 - 



2cos(coAT) 



AT 2 



This means in particular that 



m,\y{(x>r\)u?r\ = eie 2 



2 cos (<x>AT) 
AT 1 



which by (4) allows to solve for AT according to 



AT 2 = 



1 + ^ m ir(^i)^i 




(44) 



for certain values of r 2 and a>. 

We can choose e 2 £ R such that eie 2 < 0, r 2 > such that (42) holds true, and wel\(0) with 
\a>\ sufficiently small such that (44) is well-defined and < u)AT < |. This proves the claim. □ 

5 What can we learn from the approximate model? 

We emphasize that the presented construction of solutions relies sensitively on the simplicity of the 
force field F in (7), which was chosen such that it has a global inverse I as defined in Lemma 2.5. 
Any generalization of this force field that does not have a global inverse, in particular the one of 
WF, will require a new technique. However, the approximate model uses the same delay function 
(4) that is also used in WF. Therefore one can expect that many mathematical structures appearing 
in the approximate model will also arise in WF. As examples we point out that the data needed to 
define the energy functional discussed in Section 3 is exactly the same as the one needed to define 
the corresponding energy functional in WF (compare figure 3 in [WF49]), and the stable orbits 
actually coincide with the ones in WF (compare [Sch63]). Given these similarities it seems rea- 
sonable to expect that also in WF the natural choice of initial data that uniquely identifies solutions 
will be of the same type as in the approximate model given in Definition 2.1, i.e. Figure la. In this 
respect it is comforting to note that the considered initial data is already sufficient to compute the 
energy functional. The only additional information in the initial data considered here, i.e. (10) and 
(11), is about how smooth the solutions are. Smoothness will become a more delicate issue when 
considering the full WF interaction as (2) also depends on the acceleration. 

Concerning the dynamics of many charges the question of typical behavior becomes relevant 
and a generalization to many particles of the approximate model is a good candidate for studying 
measures of typicality for delay dynamics of this kind. Note that while the generalization of the 
above uniqueness result requires a slightly more sophisticated proof the results of Section 3 and 
Section 4 have straight-forward generalizations to many charges. 
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